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I.  INTRODUCTION 

It  is  a  common  knowledge  that  plates  are  essential  structural  elements 
for  constructing  marine  vehicles.  To  serve  their  functions,  most  plates 
may  have  holes  and  elastic  inclusions  and  be  subjected  to  temperature 
gradients.  It  is  the  objective  of  this  research  to  develop  an  approach  for 
the  determination  of  temperature  distributions,  thermal  stresses  and  de¬ 
formation  of  rectangular  plates  with  circular  inclusions.  It  is  assumed 
that  both  the  plates  and  inclusions  are  individually  homogeneous  and 
isotropic. 

In  the  last  decade,  several  papers  on  membrane  thermal  stresses  in 
infinite  plates  with  holes  appeared  in  the  literature.  Goodier  and  Florence 
published  a  number  of  articles  [1-6]  on  localized  thermal  stress  at  holes, 
cavities  and  inclusions  ca*’  *-ed  by  disturbed  uniform  heat  flow.  It  is  well 
known  that  in  an  infinite  plate  a  steady  heat  flow  under  a  uniform  temperature 
gradient  does  not  induce  thermal  stress  provided  the  plate  is  free  to  expand. 
Thermal  stress  is  induced,  however,  if  the  uniform  flow  of  heat  is  disturbed 
by  the  presence  of  an  insulated  circular  hole,  or  by  an  inclusion  of  another 
material.  The  influence  of  the  shape  and  rigidity  of  an  elastic  inclusion  on 
the  transverse  flexure  of  thin  plates  was  investigated  by  Goland  [7], 

Tauchert  [8]  analyzed  thermal  stresses  at  spherical  inclusions  in  uniform 
heat  flow.  Recently  Hoffman  and  Ariman  [9],  Rao,  Rao  and  Ariman  [10, 

11]  have  used  the  least  squares  boundary  point  matching  method  to  in¬ 
vestigate  the  temperature,  bending  and  membrane  stress  distributions 
for  a  rectangular  elastic  plate  with  circular  [10]  and  elliptic  f  11]  holes 
subject  to  heat  flow. 
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On  the  other  hand  in  many  cases  a  satisfactory  solution  for  a  given 
boundary  value  problem  can  be  obtained  by  the  method  of  point  matching. 

A  simple  description  of  this  technique  is  that  solutions  are  used  which 
satisfy  the  governing  differential  equations  exactly,  while  approximating 
the  boundary  conditions  at  discrete  points  on  the  boundary.  This  proce¬ 
dure  proves  very  useful  for  problems  in  which  no  exact  theoretical  analysis 
can  be  formulated.  Furthermore  a  single  computer  program  can  be  used 
for  any  problem  governed  by  that  equation  without  regarding  the  shape  of 
the  boundary. 

The  concept  of  the  least  squares  boundary  point  matching  method 
(sec  Appendix  A)  occurred  to  several  investigators  as  early  as  1934, 

Slater  [12]  and  1937,  Barta  [13].  The  method  was  not  used  extensively 
until  I960  when  Conway  [14],  after  recognizing  the  method's  potential  for 
a  wide  class  of  problems,  named  this  technique  "point-matching".  Since 
then,  the  method  has  been  developed  particularly  for  the  bending  of  plates 
by  Lcissa  and  Niedenfuhr  [15]  and  has  been  applied  to  numerous  other 
boundary  value  problems  [16-19].  Eringen  and  his  co-workers  extended 
this  method  to  stress  concentration  problems  in  shells  with  holes  and  to 
the  problem  of  two  normally  intersecting  circular  cylindrical  shells  [20-22]. 

In  the  first  problem  of  this  study,  the  temperature  distribution  in  a 
finite  rectangular  plate  with  a  circular  inclusion  is  investigated  by  the  same 
method.  The  plate  is  heated  uniformly  through  the  thickness,  that  is  the 
temperature  varies  only  along  the  coordinates  within  the  plane.  Two 
opposite  edges  of  the  plate  are  kept  at  constant  temperatures  while  the 
other  edges  are  insulated. 


In  the  second  problem,  the  ueflcction  and  moment  distribution  in  the 
same  plate  model  with  simply  supported  edges  is  investigated  by  the  least 
squares  point  matching  method.  The  heat  flow  consists  of  an  arbitrary 
temperature  variation  across  the  thickness  of  the  plate,  but  with  no  varia¬ 
tion  over  the  plate  surface.  For  both  problems  the  solutions  are  derived 
in  the  form  of  an  infinite  series  in  terms  of  circular  cylindrical  coordinates 
r  and  d.  The  boundary  conditions  at  the  inclusion  interface  are  satisfied 
exactly  while  the  ones  at  the  simply  supported  edges  of  the  rectangular 
plate  are  satisfied  at  selected  points  by  the  least  squares  method.  Numeri¬ 
cal  results  are  presented  for  temperature  distribution,  deflections  and 
moments  as  a  function  of  the  ratio  of  elastic  properties  of  the  plate  and 
inclusion,  and  the  inclusion  size. 


II.  TEMPERATURE  DISTRIBUTION 


A  rectangular  elastic  thin  plate  of  uniform  thickness  is  subjected  to 
heat  flow  in  the  x-direction  (Fig. 2.  1.  )The  edges  at  y  =  t  b,  of  the  plate 
are  insulated.  The  surfaces  at  x  =  ±  a  are  held  at  constant  temperatures 
Tj  and  T^  with  Tj>  T^.  A  circular  cylindrical  inclusion  with  different 
thermal  properties  than  the  plate  is  centrally  located  with  the  plate.  The 
governing  equation  for  the  steady-  state  heat  conduction  in  a  thin  plate  with 
no  sources  or  sinks  is  the  Laplace  equation  [23]. 


where 


V2  T  =  0 


v2  ,  .31  +  I  JL  + 

ar2  r  3r 


(2. 1) 


1 

2 


(2.  2) 


r  38 

and  T(r,  8)  is  the  temperature.  The  boundary  conditions  are  given  by 


T  =  Tj  at  x  =  a 
T  =  T2  at  x  =  -a 


(2.  3) 


(2.  4) 


ST 

dy 


=  0  at  y  =  ±  b 


STa 

T  =  T,  and  K 
a  b  a  Sr 


ST, 


=  K, 


at  r  =  R 


(2,5) 


b  Sr 

Equations  (2.  3)  and  (2.  4)  are  the  boundary  conditions  at  the  edge  of 
the  plate,  while  equations  (2.  5)  apply  at  the  boundary  of  the  circular  in¬ 
clusion  and  the  plate.  In  these  eqs.  K&  and  represent  thermal 
conductivity  of  the  inclusion  and  plate  materials  respectively. 

The  general  solution  to  Eq.  (2.  1)  in  plane  polar  coordinates  can  be 
represented  as 


1  n  -n 

T(r,  6)  =  4  (a  *i*  3  in  r)  +  £  [a  r  +  g  r  )  cos  n6 

2  o  o  n=i  n  n 


+  (y  rn  +  5  r  n)  sin  n8l 
n  n  ■* 


(2.  6) 


where  a  ,  S  ,  v  and  6  are  unknown  coefficients  to  be  determined  from 
n  Kn  n  n 

the  boundary  conditions.  By  considering  thy  symmetry  of  the  problem 
about  the  x-axis,  this  solution  reduces  to: 

For  the  plate:  Th(r,9)  =  ^  +  ^  [(anKrn  +  3rKr”n)  cos  nS  (2.7) 


n=  1 


nb 


nb 


i  00 

For  the  inclusion:  T  (r,8  )=:y  (a„_+ £  _in  r)+  £  (a  rn+ g  r_n)cos  n6  (2.  8) 


oa  oa 


n=l 


na 


na 


whe^e  T^  and  Ta  represent  temperatures  in  the  plate  and  in  the  inclusion 
respectively.  Ta  should  remain  finite  as  r  approaches  zero  in  the  in¬ 
clusion,  therefore  it  is  easily  seen  that  8  =8  =0  and  T  takes  the 

7  oa  Kna  a 


following  form 


c 


Ta(r’6)  =  2 


oa 


+  £  (a  rn)  cos  n  8 


n=l 


na 


(2.  9) 


Eqs.  \Z.  7)  and  (2.  9)  by  the  use  of  boundary  conditions  around  the  circular 
inclusion  Eqs.  (2.4)  and  (2.5),  yield  the  following  relationships: 
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“ob  +  eob  K  =  "oa 

(2.  10) 

%bRn+enbR'n=%aR" 

(2.  11) 

fob'0 

(2.  12) 

K,  (n  .  Rn_1- S  ,  E'"'1)  =  K  (a  R11'1) 
b  nb  Fnb  '  a '  na  ' 

(2.  13) 

From  Eqs.  (2.  10)  and  (2.  12) 

aob  =  aoa 

(2.  14) 

anf  from  (2.  11)  and  (2.  13) 

,  K, -K 

n  _  _  p2n  b  a 

Pnb~anbK  K  +K, 

a  b 

(2.15) 

Substituting  Eqs.  (2.  12),  (2.  14)  and  (2.  15)  into  (2.7)  and  (2.9) 

one  obtains 

a  ,  ®  ? 

for  r>R  T, (r,  0)  =  +  V  a  ,(rn+R<lnY  r~n)  cosn0 

b  2  n=  1  no 

(2.  16) 

a  ,  °° 

for  r  <  R  T  (r,  0)  J.  2J  a  .  (1  +  Y)  rn  cos  n0 

a  2  n=  j  nb 

(2.  17) 

whe  re 

K, -K 

b  a 

i?r  — 

•  ku+k 

b  a 

(2.  18) 

For  the  case  of  an  circular  hole  in  the  same  plate  K  -  0,  Y  =  1  and  Eq. 

a 

(2.  16)  is  in  agreement  with  the  one  given  by  Rao  et  al  [10], 


The  unknown  coefficients  ao^  and  a  ^  are  now  determined  by  satisfying 
the  boundary  conditions  at  the  outer  edges  of  the  plate  (Eqs.  (2.  3)  and  (2.4)). 
However,  since  the  boundaries  x  --  ±a  and  y  =  ±b  are  not  coordinate  lines 
for  the  temperature  function,  an  exact  solution  cannot  be  found.  To  obtain 
an  approximate  solution  the  infinite  series  in  Eqs.  (2.  16)  and  (2.  17)  is 
truncated  at  n  =  N  terms,  Having  n  +  1  coefficients  to  be  determined.  M 
points  are  chosen  along  the  plate  outer  boundaries.  As  there  is  one 
condition  for  each  point,  M  equations  are  generated.  For  M  =  N  +  1,  the 


3 


equations  may  be  solved  for  the  required  coefficients.  However,  if 
M  >  N  +  1  and  the  unknowns  are  determined  in  the  least  squares  sense, 
then  as  M  -*«  (his  process  becomes  that  of  minimizing  the  integral  of 
(he  squared  error.. 

In  order  to  get  some  numerical  values  for  the  problem  and 
o  ° 

are  chosen  as  150  F  and  100  F  respectively.  The  number  of  boundary 
points  taken  is  forty  (M=40)  Fig.  (2.  2)  and  eight  terms  (N=8)  are  considered 
in  the  series.  Thus  forty  equations  for  nine  unknowns  are  generated. 

Solving  these  in  the  least  squares  sense  yields  the  necessary  coefficients 

V 

In  Fig.  (2.  3),  the  variation  of  the  nondimen sionalized  temperature 
T-  T2 

function  0  =  =  —  is  presented  on  the  perimeter  of  the  inclusion  in  a 

T1'T  2 

square  plate.  The  curves  are  given  for  various  sizes  of  the  circular  in¬ 
clusion  and  no  inclusion  case  with  R/b  =  0.  The  conductivity  of  the  in¬ 
clusion  is  assumed  to  be  one-tenth  of  the  plate  \  =  .  1.  It  is  seen  that  as  the 
size  of  the  inclusion  becomes  larger  the  temperature  at  the  inclusion  inter¬ 
face  increases  and  approach  the  constant  values  at  the  outer  edges.  In 
Fig.  (2.  4)  the  inclusion  size  is  held  constant  (R/b  =  .  5)  and  the  conductivity 
ratio  X  is  varied.  Lowering  the  conductivity  ratio  increases  the  tem¬ 
perature  at  the  interface  to  the  case  of  an  insulated  circular  hole,  in¬ 
creasing  the  conductivity  ratio  X  approaches  the  case  of  the  infinitely  con¬ 
ducting  inclusion  which  holds  all  of  the  heat.  Figure  2.  5  displays  the 
effect  of  plate  size  on  the  variation  of  normalized  temperature  0  for  a 
given  inclusion  size  R/b  =  .  5  and  a  prescribed  conductivity  ratio  X  =  .  1. 

It  is  seen  that  for  a/b  =  4  temperatures  on  the  perimeter  of  the  inclusion 
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Si 

are  quite  lower  than  the  case  of  ig  =  .  25. 

Figure  2.  6  shows  the  normalized  temperature  distribution  for  a  square 

R 

plate  with  a  circular  inclusion  determined  by  ^  =  .  5  and  \  =  10.  The  curves 
represent  the  temperature  variation  for  various  values  of  y  =  As  ^ 
increases  the  temperature  approaches  the  linear  distribution,  thus  verifying 
the  fact  that  the  effect  of  the  inclusion  dies  out  as  the  distance  from  the 
hole  increases.  The  similar  temperature  distribution  is  presented  in  Fig. 

2.  7  for  a  lower  conductivity  ratio  X  =  .  1.  Again  the  temperature  variation 
is  almost  linear  for  y  =  1.  0  while  the  x-axis  (y  =  0)  shows  the  effect  of 
the  inclusion  most.  However  due  to  the  lower  conductivity  ratio,  the  tem¬ 
perature  drop  from  the  edge  of  the  plate  to  the  inclusion  is  less  than  in 
Fig.  2,  6.  Finally  for  square  plate  with  a  given  inclusion  size  of  .  75, 

Fig.  2.  8  shows  the  temperature  distribution  as  a  function  of  x,  while  the 
conductivity  ratio  is  varied.  Again  the  temperature  differential  from  the 
plate  edges  x  =  ±a  to  the  inclusion  interface  is  les3  for  the  lower  con¬ 
ductivity  ratios. 

IH.  THERMAL  BENDING 

In  the  second  part  of  this  study  thermal  bending  problem  of  a  rectangular 
elastic  thin  plate  with  an  elastic  inclusion  by  the  least  squares  point  matching 
method  is  investigated.  The  upper  and  lower  surfaces  of  the  plate  are  held 
at  constant  temperatures  T^  and  T^  respectively.  (Fig.  3.  1)  The  temperature 
is  assumed  to  vary  in  the  z-direction  only.  Hence 

T  =  T(z)  =  T  +  T  z 
o 


(3.  1) 


where 


T,'+T  ' 


T,'-T' 


T  =  with  t:  >  T'  and  t  =  —  v  ■  • 

o  2  12  2 


(3.2) 


The  governing  equation  for  the  transverse  displacement  w  of  a  plate 
of  uniform  thickness  t  subjected  to  distributed  transverse  loading  and 


nonuniform  heating  is  [  24] 


„  2  2  ,  1  „2  „ , 

D  7  V  w  =  p(x,  y)  -  -  V  Mt 


(3.  3) 


where 


p(x, y)  =  Distributed  transverse  loading  per  unit  area. 

3  2 

D  =  Bending  rigidity  of  the  plate  per  unit  length  =  E  t  /12(l-v  ) 

.  t/2 

=  Thermal  bending  moment  =  a.E  T  (z)  dz 

1  “-t/2 

v  =  Poisson’s  ratio 


by 

E  =  Young’s  modulus 

at  =  Coefficient  of  thermal  expansion 

For  the  case  of  no  transverse  loading  and  a  temperature  distribution 
described  by  Eq.  (3.  J),  the  governing  equation  reduces  to 


2  2 

V  V  w  =  0 


(3.  4) 


The  solution  of  equation  (3.  4)  in  plane  polar  coordinates  can  be  represented 
as  [25] 

2  2  - 1  3 

w  =  A  +  B  Xnr  +  C  r  +  D  r  /nr  +  [A,  r8  +  B.  r  +  C.  r 
o  o  o  o  1  1  1  1 

- 1  3 

+  D^r  /nr  ]  cos  0  +  [A'  r  0  +  Bj  r  +  Cj  r  +  Dj  r/nr]  sin0 

+  Drt  r20+  S  (A  rm+  B  r"m+  C  rm+2+  D  r'm+2)  cos  m0 
o  _  ~  ~  m  m  m  m 

m**4i  f  • « • 

+  S  <  Am  *m  +  *’m  r'm+  rm+2  +  sin  m  0  (3.  5) 

m=2,  3*  •  •  111 

By  considering  the  symmetry  of  the  problem  two  solutions,  one,  using  the 


subscript  b,  for  the  plate  and  one  using  the  subscript  a,  for  the  inclusion 
are  formulated  from  Eq.  (3.  5) 

For  the  plate: 

wb  =  Aob+  Bob  *nr  4  Cob  r2  +  Dob  r2  tor 


+  T.  (A  ,  r 
_  .  '  mb 
m=2,  4, . . . 


(A  ,  rm+B  .  r-m+C  .  rm+2+  D 
'  mb  mb  mb  r 


-m+  2.  - 

r  )  cos  m  0 


For  the  inclusion: 

w  =  A  +C  r2+  D  r2inr+r  (A  marm+  C  rm+2)cosm8 
3  03  03  03  m=2, 4, . . .  ma  (3.7) 

Here  A0a’  Aob  to  Doa’  Dob  and  Ama’  Amb  to  Dma>  Dmb  rePresent  the 
unknown  coefficients  to  be  determined  from  the  boundary  conditions.  At 

the  inclusion  interface  (r=R)  there  are  four  boundary  conditions  to  con¬ 
tend  with.  They  are 

<Mr>a  =  <Mr»b;  (Qr»a  =  <Qr>b 


w  =  w,  ; 
a  b 


(3.  8) 


where 


.  .  ^  r  a  w  ,  ,  i  aw,  i  aw, .  mt  n 

Mr  =  radial  moment  =  -D!  — j— +  V  ( -^  — J-+  f  +(T^)B  j 

*  3r  3M  r  36 

Q  =  Kirchoff  shear  force  =  Q„  +  —5-5 — 
r  r  o  0 


Q  =  vertical  shear  force 
r 

Mrg  =  twisting  moment 


Substituting  Eqs.  (3.  6)  and  (3.  7),  into  (3.  8)  one  obtains  the  following  four 
equations  by  combining  the  terms  not  under  the  summation  Riga. 
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DJZ  Coa(1+va)+  Doa(3+2^nR)+  DQa va (1+2  in  R )  + 


MTa  1  r  Bob 


(1-v  )D.  J  '  °b  ^  2  (V1,+  2  C„b(1+V  + 

n 


M, 


Dob  (3t2  *"R>  +  nob  V1+21nR>  I  -(]  ffi 

b  b ' 


4D 


D 


oa 


4  D 


a  R 


=  D, 


ob 


(3.  9) 
(3.  10) 


Aoa+  Coa  r2+  Doa  r2  *nR  =  Aob+  Bob  *nR  +  CobR‘ 
+  DQb  R2  inR 


(3.11) 


2RCoa+Doa(R+2R*nR>  = 


B 


ob 


-r—  +  2RCob+  Dob(R+  2R  *nR)  (3‘ 12> 


These  four  equations  may  be  solved  for  A  ,  ,  B  ,  ,  C  ,  and  D  . 

ob  ob  ob  ob 


D  ,  =  n  D 
ob  oa 


Bob  =  al  Coa+  a2  Doa+  a3 


Cob  a4  <2oa+  05  Doa+  a6 


ob 


=  a?  Aoa+a8  Coa+a9  °oa+ a10 


(3.  13) 


where 


R  r  <  i+vb)  -  tki+ve)] 


a- 


a 


aJ(l+2  inR)/2] 
-R2/ 2  r  MTa 


M 


Tb 


d-Va)Da 


(1_Vb,Db  J 

a4  =  fn(l  +  va)+  (1  -  vb)]/2 

a5  =  f(l+2inR)/4]  [  (l-vb)  -  n (l-va) ] 


a. 


=  a3/-2R' 


a?  =  1.0 


a 


8 


=  [R2  -  inR  ax  -  R2  a4] 


aQ  =  fR2  inR  (1-t)  -  inR  a,  -  R2  a.l 


IT  1  A' 


Combining  the  terms  under  the  summation  sign  in  the  substitution  of 
Eqs.  (3.  6)  and  (3.  7)  into  {3.  8)  the  following  four  equations  in  six  unknowns 


Then  Eq.  (3.  6)  can  be  written  as 


w,=  F.A  +  F-C  +  F-.D  +E„+r  (§_A  +  §,C  )cosm0 

b  M  oa  *2  oa  -3  oa  -4  ,  .  ...  5  ma  ’6  ma' 

IT*  **Li|  *  f 

(3.  18) 


where: 


§j-  =  a?,  ?2  =  ag  +  aj  in  r  +  a4  r 

2  2 

E,  =  a^  +  a2  in  r  +  a5r  +n  in  r 


?4  =  a10+a3  *nr  +  <V 


m 


m  .  „  m+2 


-m+2 


F5  =  0?  r  +  g5  r  *“+  gj  r  +  03  r 

_  .  m  ,  _  -m  ,  „  m+2  .  -  -m+2 

?6  =  e8  r  +  86  r  +  8,  r  +  r 


(3.  19) 


'2  ‘  r  »4 

The  condition  that  the  outer  edges  of  the  plate  be  simply  supported  requires 
that 


at  x  =  ±  a  and  y  =  ±  b: 

a2w  S2w,  M 

wb  =  0  and  Mn  =  -Db(— £-+  *b  ~£T >  ’  *  0  <3' 20) 

where  n  and  t  denote  normal  and  tangential  directions  at  the  plate  boundary 
» <_-.->pectively. 

Substitution  of  Eq.  (3.  18)  into  the  outer  edge  boundary  conditions  re¬ 
sults  in  an  infinite  number  of  equations  for  the  infinite  number  of  unknown 
coefficients.  As  in  the  case  of  the  first  problem,  the  least  squares  boundary 
point  matching  method  was  employed.  Namely  the  simply  supported  boundary 
conditions  are  satisfied  at  selected  points,  along  the  plate  outer  boundaries. 

In  order  to  present  som  e  numerical  results,  Tj  and  T^  are  chosen  as 
150°F  and  100°F  respectively  while  the  plate  thickness  is  Forty 

points  are  taken  along  the  outer  boundaries  of  the  plate  (Fig.  2.  2)  and 
evght  terms  are  considered  in  the  series.  Thus  eighty  equations  in  eleven 
unknowns  are  generated.  Solving  these  in  the  least  squares  sense  yields 


turn 
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the  necessary  coefficients. 

w 

In  figures  3.  2  and  3.  3  the  deflection  function  — g-  is  presented  as  a 

a 

function  of  X  =  —  .  The  ratio  of  the  inclusion  rigidity  to  plate  rigidity, 

D  a 

r  =  -jr-  ,  is  chosen  as  .  1  and  10  respectively  and  the  prescribed  value  of 

b  MTa 

the  ratio  of  the  thermal  moments  p  =  — —  =,  is  1.0.  Each  figure 

MTb 

presents  the  variation  of  the  deflection  function  for  three  different  inclusion 
sizes  -^  =  .  25,  .  50,  and  .75.  In  Fig.  3.  2  the  inclusion  is  less  rigid  than 
the  plate  and  it  is  seen  that  the  larger  inclusion  increases  the  plate 
deflection.  However,  in  Fig.  3.  3,  the  inclusion  is  more  rigid  than  the  plate 
and  increasing  the  inclusion  size  decreases  the  plate  deflection. 

In  Figs.  (3.  4)  and  (3.  5)  the  deflection  function  is  presented  on  the 
perimeter  of  the  inclusion.  The  rigidity,  thermal  moment  and  Poisson 
ratio  and  inclusion  size  are  held  constant  while  the  deflection  at  various 
radii  is  plotted.  A  comparison  of  these  two  figures  indicates  that  the 
increase  in  the  rigidity  of  the  inclusion,  cause-*  lower  deflection  values  at 
corresponding  radii. 

A  summary  of  the  effect  of  the  inclusion  rigidity  is  presented  in  Fig. 

3.  6.  In  this  figure  for  a  prescribed  inclusion  size  the  deflection  function 
variation  along  the  x-axis  is  shown  as  a  function  of  the  inclusion  plate 
rigidity  ratio  -p.  The  rigidity  may  be  increased  to  the  limit  of  the  rigid 
inclusion  where  the  deflection  is  constant  across  the  inclusion  surface. 

At  the  other  extreme,  the  rigidity  is  lowered  to  the  point  where  the  inclusion 
gives  no  support  and  acts  as  a  circular  hole. 

Finally  in  Fig.  3.  7  the  effect  of  the  ratio  of  thermal  moment  p  on  the 
variation  of  the  plc.te  deflection  is  displayed.  For  a  given  inclusion  size, 
constant  elastic  properties  (Ea,  E^,  and  v^),  the  increase  in  the  ratio 


MNH-  7.'l-2 


-  14  - 


P  winch  is  directly  related  to  the  ratio  of  thermal  expansion  coefficients 

a.t 

•—  ,  causes  larucr  values  in  the  deflection  function. 

"b 


uni*-  /:?-:> 
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APPENDIX  A 

THE  LEAST  SQUARE  BOUNDARY  POINT  MATCHING  METHOD 


Given  a  linear  function  of  the  variables  g.  defined  by 
N 

(A-l) 

i=l 

where  the  g.  may  be  functions  of  any  finite  number  of  coordinates  (x^). 

It  is  desired  to  determine  the  coefficients  a.  such  that  f(x.)  satisfies 
some  specified  conditions  on  some  contour  L.  If  L  is  not  a  coordinate 
line  of  the  coordinate  system  in  which  the  x^  are  defined,  it  is  presently 
impossible  to  do  this  exactly.  However,  if  at  some  finite  number  M  of 
points  on  L  the  specified  conditions  are  satisfied,  the  following  system 
of  equations  results. 


-  N 


(A- 2) 


Fm  L«*k’j  '  FmLs?,  ai«i  <V  j  *  bm 

where  (m  =  1,...,M).  If  F  is  a  linear  operator 
N 

b  =  E  a  F  (g  (x.))  (m  -1 . M)  (A- 3) 

m  j_  j  i  m  i  k 

Then  if  M  =  N  the  system  is  determined.  However  if  M>  N  consider  a 
solution  to 

M  „  N  2 


fbm  -  E,  *iFm(gi  (*k»  »  j  =d 
m=l  i=l 


(A- 4) 


where  d  is  minimized  with  respect  to  the  desired  solution  vector  a.. 


That  is 


rbm-r  aiFm,g.,xk)sr 

a.  m=l  i=l 

i 


It  can  be  shown  that  if 


a.  =  a 
i  — 


Fm'gi(xk»>=  tCl 


iirm-  n-2 


1.4  - 


and  b  =  b  then  for  a  to  minimize  {A- 2)  for  M  >  N  it  is  necessary  that 
a  be  a  solution  of 

[C]TfC]  a  =  [C)Tb 
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Fig,  3,1  Simply  supported  rectangular  plate  with  circular  inclusion 
coordinate  system. 


deflection  as  a  function  of  0  wit! 


